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Abstract. The possibility of a connection between dark energy and gravity through 
a direct coupling in the Lagrangian of the underlying theory has acquired an increasing 
interest due to the recently discovered capability of the extended quintessence model 
to encompass the fine-tuning problem of the cosmological constant. The gravity 
induced "R-boost" mechanism is indeed responsible for an early, enhanced scalar field 
dynamics, by virtue of which the residual imprint of a wide set of initial field values is 
cancelled out. The initial conditions problem is particularly relevant, as the most recent 
observations indicate that the Dark Energy equation of state approaches, at the present 
time, the cosmological constant value, wde = — 1; if confirmed, such observational 
evidence would cancel the advantage of a standard, minimally coupled scalar field as 
a dark energy candidate instead of the cosmological constant, because of the huge fine 
tuning it would require. We give here a general classification of the scalar-tensor gravity 
theories admitting R-boost solutions scaling as a power of the cosmological redshift, 
outlining those behaving as an attractor for the quintessence field. In particular, 
we show that all the R-boost solutions with the dark energy density scaling as the 
relativistic matter or shallower represent attractors. This analysis is exhaustive as for 
the classification of the couplings which admit R-boost and the subsequent enlargement 
of the basin of attraction enclosing the initial scalar field values. 
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1. Introduction 

In the last few years, cosmology has gone through a deep revision of the basic 
ideas on which it was used to rely. In particular, the most puzzling problem 
for the old, "standard" CDM scenario, arised when the observations of distant 
la Supernovae evidenced an accelerated expansion of the Universe, through the 
reconstructed magnitude- redshift relations [TJ |2] . Since these very early indications, the 
case for a "Dark Energy" component able to accelerate the cosmic expansion became 
increasingly stronger when the observations of the CMB anisotropies, together with large 
scale structure data, clearly revealed a very close-to- flat, low-density Universe [Hill]. 
Many phenomenological explanations for the dark energy involve classical, minimally- 
coupled scalar fields violating the strong energy condition ("quintessence" fields, 0-[E] 
and references therein), or a phantom energy (see, e.g., [Oj) which, as a cosmological 
constant, violates the weak energy condition. Most of these models, however, suffer the 
worrying problem of initial conditions, in the sense that an accurate and sometimes 
unphysical tuning of the initial values of a field is required in order to reproduce 
the cosmological conditions observed today (in particular, the dark energy equation 
of state wde and the density parameter of the dark energy component). This problem, 
affecting almost all dark energy models involving scalar fields, is partially alleviated 
when, depending on the potential on which the field is assumed to evolve, the scalar 
field equation of motion admits tracking solutions: in that case, indeed, the present 
value of the cosmological parameters, as requested by the observations, can be reached 
starting from some more or less extended range of initial values of the field. However, 
even for those models admitting a tracking behavior, the problem of initial values is 
now becoming particularly serious, because the observation bounds on the dark energy 
equation of state are increasingly converging towards a value of wde very close to —1. It 
is known that, the closer wde is to the cosmological constant value, the smaller the range 
of allowed initial value for the corresponding scalar field has to be, since the dynamics 
of such a field is extremely constrained by the flatness of the potential in which the 
field evolves [TU]. Furthermore, the best fit of the latest Sn la data (the Gold dataset, 
[TT] ) includes values w^e < — 1? which cannot be obtained in the context of minimally- 
coupled quintessence models. In light of these considerations, it is worth to point our 
attention to extensions of General Relativity where, standing the presence of a scalar 
field as a dark energy candidate, the dynamics of the field itself can be, at early times, 
strongly modified by gravitational effects. In particular, we will focus on scalar-tensor 
theories (see, e.g., 12J where a scalar field, non-minimally coupled to Gravity, acts as 
a dark energy component at recent times A dark energy component, as it 

arises in scalar-tensor theories, has been proven to have several advantages with respect 
to minimally-coupled scalar fields: the direct coupling of the field to the Ricci scalar 
in the Lagrangian of the model, makes the field undergoing an enhanced dynamics at 
early times, known as "R-boost" , with the appealing consequence that the characteristic 
thickness of the basin of attraction for tracker solutions is preserved even for wde close to 
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— 1 [27J. Furthermore, as noticed in [22] [22] [22] [21], such an "extended quintessence" 
(EQ) component can cross the cosmological constant value, getting wbe < — 1- In 
the past however, those interesting results were obtained considering particular cases of 
scalar-tensor theories of gravity, while a general treatment is still missing. The purpose 
of the present paper is to fill this gap. We provide here a general classification of all 
the possible scalar-tensor models of Dark Energy which admit R-boost trajectories for 
which the quintessence energy density scales as a power law of the cosmological scale 
factor. In the following, we refer to those as scaling solutions. Although in a different 
context, for purpose and methodology our paper follows the footsteps of j2S] performed 
in the context of minimally coupled quintessence models. 

The plan of the paper is the following: in section [21 we describe the working framework, 
giving the basic definitions cosmological equations; in section El we search for generic 
scaling solutions of the Klein-Gordon equation, and classify the coupling functions 
according to the resulting type of scaling; in section 0] we check these solutions for 
stability, verifying which ones among the theories we found produce stable attractors; 
finally, in section we draw our conclusions. 



2. Framework in generalized theories of gravity 



S = d^xyf^g 



- f(<f>, R) - -u;( - V{<f>) + £ fluid 



(1) 



The class of Generalized Theories of Gravity we refer to is described by the action 

where g is the determinant of the background metric, R is the Ricci scalar and <fi 
is a scalar field whose kinetic energy and potential are specified by u((j>) and V(4>), 
respectively. £fl U id includes contributions from all components different from and 
K = 8nG* represents the bare gravitational constant. We also work in natural units, 
c = 1. The action above has been first considered in a cosmological context in 
[HE]- In particular, the classes of theories in which /(</>, R) assumes the simple form 
f((j), R)/2k = F((j))R/2 have been considered in the context of dark energy cosmologies 
(extended quintessence, see [22] and references therein). Note that the gravitational 
constant differs from the one measured in Cavendish like experiments by corrections 
being negligible in the limit ujjbd 1 [20] where ujjbd is defined as 



F 

UJBD = TT2 ( 2 ) 

and F^ is the derivative of F with respect to <fi. For a review of the allowed values of the 
ujjbd parameter and other constraints on generalized theories of gravity we send back 

to [27]- [ID]. 

Compared to general relativity, the Lagrangian has been generalized by introducing 
an explicit coupling between the Ricci scalar and the scalar field, achieved by replacing 
the usual gravity term R/lQirG with the function R)/2k. This new term, which 
has the effect of introducing a spacetime dependent gravitational constant, may either 
be interpreted as an explicit coupling between the quintessence field and gravity, or as 
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a pure geometrical modification of general relativity admitting a non-linear dependence 
on R. 

In the assumption of flat cosmologies, the line element can be written as ds 2 = 
— dt 2 + a 2 {t)b~ i jdx % dx : > , where a(t) is the scale factor and t represents the cosmic time 
variable; the expansion of the Universe and the dynamics of the field are determined by 
Friedmann and Klein Gordon equations 

H 2 = (-) 2 = i ( P fluid + U 2 +V(cf>)- 3HF) , (3) 



a) 3F 

+ 3tf = - V , (4) 

where dot means derivative with respect to the cosmic time t, and V§ are the 
derivatives of the coupling F and of the quintessence potential V with respect to 0. 
In this context, the Ricci scalar can be written in terms of the cosmo logical content of 
the Universe 



-4 



P fluid + %P fluid + r-W + 3F + 9HF , (5) 



where pfi U id and pfi U id are the energy density and pressure summed up over all possible 
cosmological components but 0. 

As it was shown in j^Sl, all species but satisfy the usual conservation equations 
pi = —3H(pi + Pi), where pi and Pi are respectively the energy density and pressure 
of the i-th component. Unlike minimal coupled models, in extended quintessence the 
energy density 

P<t> = ~0 2 + V{4>) (6) 

is not conserved, i.e. it does not obey the relation p ( f,-\-3H(p < f > -\-p ( j > ) = 0; nevertheless, one 
can still define a conserved expression for the energy density including the contributions 
from the non-minimal coupling which reads as 



p cons. = 1^2 + _ 3^ + 3^2 (I - f) . (7) 

The latter expression can be generally very different from ©, mostly because of the 
last term, proportional to the cosmological critical density and active whenever the 
theory differs from general relativity. Although that difference is small to match 
the observational constraints, the presence of H 2 makes it relevant. This may have 
important consequences for the dynamics of the dark energy density perturbations, 
leading to effects like gravitational dragging [13]. However, as we already stressed, 
in the present work we look for trajectories of the scalar field as a function of time, 
determined by the effective gravitational potential in the Klein Gordon equation. In the 
energy density, this contribution is entirely kinetic [18] . which therefore is the relevant 
energy density component for our purpose, as we see in detail in the next Section. In 
other words, we are interested in attractor paths which solve the Klein Gordon equation 
for the field, which is unique regardless of the definition of the energy density. In this 
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perspective, it is convenient to define a third expression related to the dark energy 
density [TH] : 

P^ = \^ + V{<P)-?,HF . (8) 

The quantity above has a direct intuitive meaning for the dynamics of cosmological 
expansion, as it regroups the scalar field terms which compete with pf[ U id in the 
Friedmann equation (JSJ). 

The cases in which F{4>) has a quadratic form (both induced gravity and non 
minimal coupling theories [XT] ITS] ) or an exponential form [3U] have been widely 
investigated and reached their main achievement in what concerns the early dynamics 
of the dark energy field [21] • It was shown that as soon as there are non-relativistic 
species in the radiation dominated era, the new term containing R in the Klein Gordon 
equation has a diverging behavior 

, (9) 
where p m nro is the present value of the energy density of the components which are 
non-relativistic at the time in which the R-boost occurs, representing the only non 
zero contribution of the term —pfi U id + 3p fluid in eq.©; the remaining terms in © 
are negligible for a sufficiently small. The R-boost is caused by the effective, time 
dependent potential arising in the Klein Gordon equation because of the presence of 
the non-minimal coupling. The main point to be stressed in view of the following 
analysis is that, at least for the particular choices of coupling investigated up to now, 
the R-boost guaranties an attractor behavior. This is a crucial aspect for extended 
quintessence because it provides a valid alternative to minimal coupled tracker fields, 
achieving attractors by means of their standard potential j3S] [UJ; as we stressed above, 
this property may disappear if the present dark energy equation of state gets close to -1 
[TU] . On the other hand, the R-boost remains a viable mechanism to keep a large basin 
of attraction if approaches —1 [2Zj. These scenario have been recently constrained 
with the available cosmic microwave background and large scale structure data [3*8] . 

In this work we provide the general analysis required to classify the scalar tensor 
theories of gravity which have attractor solutions to the Klein Gordon equation because 
of the non-minimal coupling with the Ricci scalar. With this purpose in mind, we will 
operate in the framework described by equations (JJJ) — © but we will not fix a specific 
expression for the coupling F(<ft). On the contrary, we will try to classify the possible 
forms of the coupling which give rise to R-boost trajectories behaving as attractors in 
the early Universe. In order to proceed in this direction, we will follow closely what has 
been done in [HE] in the case of minimal coupling, where the authors classify the allowed 
forms of the true potential V(<p) in order to have scaling solutions. 



3. Classification of the couplings 



Let's consider a cosmological model described by (0) — (J5J) containing both a scalar 
field (f) coupled to gravity through the function F(<f)) and a contribution of perfect 
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fluid components with pressure and energy density generically indicated with pfiuid and 
P fluid- Assuming that in the Friedmann equation the terms involving the scalar field are 
negligible with respect to the fluid energy density 

Pfiuid > P<i> , (10) 

which is true in the epoch of early Universe we are interested in, and that the fluid 
scales as 

Pfiuid oc a~ m , (11) 

we want to find the forms of the coupling F(<f>) which admit solutions of the equation 
(j3J) for which 

P^oca"" , (12) 

where m and n are two integer positive defined numbers. 

Assuming that the eq.® for the Ricci scalar is still valid and since the potential 
V(<j)) in eq.pjl has no relevant effect up to recent times, we can rewrite the Klein Gordon 
equation in the following way: 

+ 3 #0 = I^^ . (13) 

The Friedmann equation (jSJ), assuming that the perfect fluid is the dominant 
component, allows us to calculate the behavior with time of the scale factor a(t) as 
/ + \ 2 / m 

a = ( T ) , (14) 



where i* stands for a fixed reference time in the RDE and we neglected any initial 
condition, assuming to work with t large enough with respect to t*. In eq. (|14j) we have 
assumed that the time dependence of F in the Friedmann equation is modest enough 
not to affect significantly the dependence on the scale factor of the right hand side in 
the Friedmann equation Q. Also, the last two terms in the Friedmann equation must 
be small in order to satisfy the condition (jlOj) . However, we will still verify a posteriori 
that these assumptions are plausible. Substituting in eq.(j!3j) we get: 

9= T0+-/WO ( — I -TT ■ ( 15 ) 

mt 2 V t J b 
At this stage the energy density of the field <fr is mainly given by its kinetic contribution, 
acquired through slow rolling onto the effective gravitational potential in the Klein 
Gordon equation. Thus p^ ~ 2 /2 and assuming the desired scaling behavior ([12)1 we 
obtain <ft time dependence 

oc r n/m . (16) 

We will now proceed by distinguishing the two cases m = n or m ^ n. 
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3.1. Case m = n 

If m = n then oc t -1 and integrating this expression we get 

= Aln-^ + 0, (17) 

where A is a constant with the dimensions of a field, i.e. proportional to the Planck 
mass m P = l/y^G) in our units. Substituting (JT7J) and its first and second derivatives 
in eq. (jl5|) we obtain the following expression 

r PmnrO U XTU J 

If m 7^ 3, the condition on the coupling, obtained by integrating eq.([T%]). is 

F = F jB(eC«-*.)-l)} (ig) 

where 

A 2 1 6-m ^ _ 2 / 3 

PmnroU 2 3-171 A U 

and is the value the coupling has at £*. Note that the combination p mnr otl has a 
direct interpretation in terms of the abundance of the non-relativistic components at 
the U time; indeed /Wo** oc (aja ) 3 p mnr j H 2 . 

If m = 3, which is the case of a Universe dominated by ordinary matter, eq.(JTHJ) 
becomes 

F* 2A 1 



B = —^^ , C=±(±-l) (20) 



F PmnrO ** 

and the form of the allowed coupling is 

P — P g[ PmnrO tj ^ ^ _ 



(21) 



(22) 

As expected, this is right the form of the coupling chosen in [HO] if we define A = 
Pmnro ^ ^gj-g ^ j g an adimensional constant and nip = 1/y/G is the Planck mass in 
natural units. In the regime <fi 3> 0*, the coupling becomes exactly the one exploited in 

[SO]. 

Notice that in minimal coupling theories the case m = n corresponds to a scenario 
in which the dark energy scales as the dominant component, thus never achieving 
acceleration unless such regime is broken by some physical mechanism, as in the case 
of quintessence with exponential potential [J21 ESj • On the other hand, in scalar-tensor 
theories of gravity the case m = n is fully exploitable in its context, since it has been 
obtained precisely with the assumption (jlUj) and neglecting the true potential V in the 
Klein Gordon equation (j3J). Actually, its relevance is on the capability to provide an 
attractor mechanism when the dark energy is sub-dominant, independently on the form 
of the potential energy driving acceleration today. Indeed, it does not exclude a different 
behavior of the <fi field at present time, when the potential V starts to have a dominant 
effect on the dynamics of 0. 
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3.2. Case m^n 

Integrating eq. (|16jl in the case m 7^ n we get: 



Ait 1 '™ -U 



+ 4>* 



(23) 



where A has the dimensions of a time derivative of a field. Substituting this expression in 
eq. ()15j) . together with its first and second derivatives, we obtain the following condition 

\ B 



where 



F 



C 




1— 



2A 2 1 (m-n)(6-n) 



5 



6 — m — n 



m — n 



(24) 



(25) 



PmnrO ^* ^ 

As in Section EIU note the combination p mn rot*^ m oc (a*/ao) 3 p mnr */ H*^ m . Again we 
have to distinguish the two cases in which the exponent B is equal or not to —1. 
For B = — 1 the condition ()24|) gives, after integration, polynomials of 

c 



1 + 



.4 



-t 



IL-l 



(26) 



Notice that if the solution (J23)) is exactly a power law of the time t, namely At* m , the 
coupling (j2Hj) becomes exactly a power law as well, as in induced gravity models (see 
e.g. JE] and references therein). Therefore, this class of gravity theories admit scaling 
solutions: however, those may arise in the radiation dominated era only, as B = —1 
with m = 4 induces n = 3; on the other hand, the same condition with m = 3 is not 
satisfied regardless of the value of n. 

For B 7^ — 1 the result of the integration of eq.(J2U) is 



l+B 



-+u- 



S + l 



- t« 



(27) 

Notice that this includes the exponential case when B = and thus m + n = 6 with 
m 7^ n. As a consequence, if m = 4 and n = 2we obtain again the coupling investigated 
in [20] for the case of radiation dominated era. The latter case may be obtained exactly 
in the limit in which t* and 0* may be neglected in (}23|) . Note that expression (|27jl is 
generalized to values of B 7^ and includes exponential functions with general coupling 
constant A contained in C and defined in (|2*3*|) . 

The same case, m = 4 = 2n yielding B = 0, also corresponds to the R-boost solution 
exploited in jTH] i n the radiation dominated era. The reason why our formalism does 
not show that the form of the non-minimal coupling in that case, F = 1/2k + £0 2 , 
is compatible with such a solution, is the following. For small values of the coupling 
constant £, the field dynamics is correspondingly reduced. Eventually one enters the 
regime in which the second term in the right hand side of ()23)1 dominates over the first 
one, yielding a scalar field value as a function of time which is effectively constant; 
this is a clearly transient phase, as eventually such regime is broken, but before that 
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the relation (|24j) remains approximately true. Although this is formally not a scaling 
solution as that found in the exponential case [3U], the values of £ may be chosen so 
small that for interesting initial conditions on 0, its variation due to the R-boost is 
not relevant at all relevant epochs in the radiation dominated era, keeping the solution 
m = 4 = In effectively valid also for the coupling considered in ^H]- We may refer 
to the solutions of the Friedmann and Klein Gordon equations in these scenarios as 
transient scaling solutions, meaning that they hold until the true R-boost dynamics 
takes over moving the field value away from its initial condition. Note finally that the 
same reasoning applies in general in all cases where F is the sum of a constant plus a 
positive power of the field, normalized by a coupling constant. 

3.3. Summary and consistency criteria 

We have found, up to now, all possible choices of the coupling F which can have scaling 
solutions verifying - ()12j) . in the sense that there are no other forms of the coupling 
admitting scaling behavior in scalar-tensor theories of gravity; for some of them, like in 
the non- minimal coupling considered in ^H] , transient scaling solutions are admitted in 
the time interval preceeding the epoch in which the R-boost motions moves the field 
away from its initial condition. We have found, as expected, exponential forms of the 
coupling [SO] which have, however, been generalized to other choices of m, n and the 
B exponent in (|27|): also, we have seen that (|26j) allows polynomials of 0. Eq. (jl9j) also 
suggests that a new family for the coupling F might be allowed, in the case m = n, 
namely made by exponentials of exponentials. Though no other coupling can allow for 
scaling solutions, we have by now no guarantee that the general solution of all these 
expressions for the coupling will indeed have an attractor behavior and this is what we 
are going to investigate in the next section. 

Before moving to that, let us briefly check the consistency of the scaling solutions we 
found with the assumptions we made. They are essentially three, concerning equations 
(JHJ), © and (JHJ). It is important to note that the variation of F induces corrections 
which are small in the limit uj jbd 3> 1. For example, F/F = (pF^/F oc 0/ ^ujjbd', 
therefore, even if the field dynamics may be important as in (|16p. the coupling constant 
may be chosen small in order to yield a small variation in time of F. More precisely, it 
is easy to see that the kinetic contribution in Q is the lowest order term in 1/lujbd] 
all the others, involving a change in F, yield terms like F^cj), which are of higher order 
due to the presence of F^, which brings another 1/ 'sqrtujBD- Indeed, it is important to 
stress again that the R-boost dynamics is caused primarily by a non-zero Ricci scalar, 
diverging in the early universe if at least one non-relativistic species is present, and 
not only to the underlying scalar-tensor gravity theory. Thus in the limit of a small 
coupling, all the three approximations mentioned above are satisfied. However, it is 
interesting to push the analysis a little further here, by computing the scaling of the 
terms in (|3|5j) coming from the scalar-tensor coupling. Concerning the last term in (JHJ), 
which may be written as —HF^/F, using (j!5|) and (j!6|) it may be easily verified that 
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it scales as 1/t m ; when compared with the first term in the right hand side of (JBJ), 
it yields the condition n < 3 both in matter and radiation dominated eras. Note that if 
one ignores the issue related with the coupling strength mentioned above, this relation 
is quite stringent, confining all the scaling solution in scalar-tensor theories of gravity 
to possess a shape not steeper than a~ 3 . The same criteria lead to n < 3 and n < 2 in 
the matter and radiation dominated eras, respectively, for neglecting the last two terms 
in the right hand side of (jSJ). Also in this case, this requirement may be bypassed by 
working in a small scalar-tensor coupling regime. 

We now turn to study which scaling solutions represent attractors for the field 
dynamics. 

4. Attractors 

Our purpose here is to check whether the various forms of the coupling F found in the 
previous section lead to attractor solutions to the Klein Gordon equation (fT3~j). With 
this aim, we will investigate whether the particular solutions found for the cases m = n 
and m n are indeed attractors or not. Following the criteria developed in [35J, we will 
proceed by linearizing the Klein Gordon equation with small exponential perturbations 
e Ar around the critical point, represented, as we will se, by our particular solution. At 
a linear level, the attractor behavior will be guaranteed whenever the perturbation will 
converge to zero with time. We shall also investigate numerically a few cases without 
linearization. 

With the following change of variable 

u = T ( 28 ) 

<Pe 

where e is the exact particular solution of eq.(fT5j). eq.(JT3J) can be rewritten as 



We will now distinguish what happens in the two cases discussed in the previous section, 
for m n and m = n. 

4-1. Attractor behavior for m ^ n 

As we have discussed in the previous section, in the case m ^ n the coupling needs to 
satisfy eq.(j24j) and the exact solution <p e depends on time as in eq.()23|). For simplicity 
in the following we will consider t 3> i* in such a way that the £* term in eq.(|23j) can 
be neglected; also, for t large enough, the initial condition 0* can be neglected too and 
eq.(|23|) reduces to = At 1_ m ; where A is an arbitrary constant as defined in eq. (|23|) . 
Substituting the first and second derivatives of 4> e in eq. (fT3*|) we get 



U(f) e + U 20 e + 3H(f) e 




2a 3 F 



(29) 



PmnrO F<$> = (6 - 7l)(m - w) e 

2a 3 F ~ m? t 2 



(30) 
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Using this expression in eq.(|29|) and with the change of variable 



we obtain 



r = lnf , (31) 

„ m-2n + 6 , (6-w)(m-n) f , , . 

u H u = , (32) 

where the derivatives are calculated with respect to r. As we can see, this equation 
admits a critical point for u — 1 and it' = 0. If we consider a perturbation 5u around 
the critical point, such that u = 1 + Su we can rewrite eq.([32|) as 

. .. m - 2n + 6 . . (6 — ra) (m — n) . /ooN 
&t" H 5u' = ^ -£u (33) 

which is a homogeneous differential equation of the second order in r with constant 
coefficients. The generic expression of the perturbation 5u is thus a combination of 
exponential terms (5u = C\e XlT + C2e A2T ) where c\ and C2 are arbitrary constants and 
A 1,2 are equal to 

m-2n + 6 |m-6| /x 

A '>' = — —±-1^ ■ < 34 > 

It is easy to see that the eigenvalues are real and negative for n <min{6,m}. As a 
consequence, for these values of the parameters, the perturbation will go to zero with 
time and the solution of the Klein Gordon equation (}3*2*|) will converge to u — 1 making 
= e to behave as a stable attractor. Notice that this discussion includes the case of 
the exponential coupling investigated in |371j . 

The numerical solution of equation (}3*2*|) . in full generality and without linearization, 
is shown in fig.Q for a somewhat large set of initial conditions. On the left side the 
time dependences of u and vl are shown, in the case (m = 4 , n = 2). Different curves 
correspond to different initial conditions for u and u' and they all converge to the stable 
attractor solution (u — 1 , u' — 0). On the right side it is shown the v! vs u plot for the 
same choices of initial conditions of the left hand side plot. 

4-2. Attractor behavior for m = n 

In the case m = n the exact solution e depends on time as in eq. ()17j) thus we have 

PmnrO _ A f 6 

2a 3 ~F ~ tf\m 

where A is the same as in expression ()17j) . With the change of variable (}3*Tj) the Klein 
Gordon equation becomes§: 

.. . f 2A 6 — m\ A6 — m.^ 

u" + u'[ — + =— (1-u) . (36) 

\(p e m J (p e m 

§ While writing the present paper, the authors recognized that eq.(47) of Pettorino et al jSOl was 
obtained under wrong assumptions. Namely, the value of <j> e found in the case m ^ n was used in 
matter dominated era (MDE) too, leading to wrong eigenvalues. 



1 , (35) 
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250 300 



T U 

Figure 1. Numerical solutions of cci. (|32[) are shown: on the left side u (black) and 
u' (red) are plotted versus r, in the case (m = 4 , n = 2) and converge to the stable 
attractor solution (u = 1 , u' = 0). Different curves correspond to different initial 
conditions for u and u': {u ,u' ) = (100,200) (solid), (100, iOO) (dashed), (250,-100) 
(dotted). The r scale is logaritmic. On the right side u' vs u is plotted for the same 
initial conditions. 



Note that, unlike what happens in the previous case, the coefficients in eq.(J3BJ) depend 
on time through <fi e {t). Nevertheless, (it = 1 , u' — 0) still behaves as a critical point 
for the equation and we are still allowed to consider a generic perturbation 5u to the 
critical point such that u — 1 + Su. Eq. (JSEI) then becomes: 

8u" + 5u — — I = — ■ 5u . (37) 

\(Pe m J <p e m 

However, we are now dealing with a homogeneous differential equation of second order 
in which the coefficients vary with time: 

5u" + P(t)5u + Q(t)5u = , (38) 

where 

, , 2A 6 — m i \ 

P{T ^M^. + — < 39) 

~, s A 6 — m , , 

«w = ^-— • (40) 

In order to find the expression of the perturbation 5u we consider the following change 
of variables: 

Su = 5ze-^I p(T)dT (41) 
in terms of which eq.(j37j) can be rewritten (if 5u ^ 0) as: 



5z" + 



Q(r) ~ \p\t) - \p\r) 



5z = . (42) 



Scaling solutions in scalar-tensor cosmologies 



13 



Figure 2. Numerical solutions of cci. (|36[l are shown: on the left side u (black) and 
v! (red) are plotted versus r, in the case (m = 3 , n — 3) and converge to the 
stable attractor solution (u = 1 , u' = 0). Different curves correspond to different 
initial conditions for u and u 1 : (uo,u' ) = (9,70) (solid), (100,20) (dashed), (50,-50) 
(dotted), (200,40) (dot-dashed). We have chosen the test values A = Mpi and 
4>o = Mpi. The r scale is logaritmic. On the right side v! vs u is plotted for the 
same initial conditions. 



It is easy to check that for our values of P and Q we get 

5z" - ^^Sz = (43) 
whose generic solution is 

(6-m) 2 

5z = ci + c 2 e~^~ T , (44) 

where C\ and c 2 are arbitrary constants. Substituting this expression in f)41j) we get the 
form of the perturbation Su: 

Clk (6-m) cok 

where A; is a constant and a = 3( * 2 ~™^2~ m * > ■ It's then immediate to see that the generic 
perturbation goes to zero for 2 < m < 6, a range including both radiation (m = 4) and 
matter (m = 3) dominated backgrounds, thus making <ft e a stable attractor. 

The numerical solution of equation (|36|) is shown in fig.© for a somewhat large 
set of initial conditions. On the left side the time dependences of u and v! are shown, 
in the case (to = 3 , n = 3) and for the test values A = Mpi and <po = Mpf. different 
curves correspond to different initial conditions for u and vl and they all converge to 
the stable attractor solution (u — 1 , u' — 0). On the right side it is shown the vl vs u 
plot for the same choices of initial conditions of the left hand side plot. 
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5. Conclusions 

The analysis of the most recent data from different cosmological probes such as Type la 
Supernovae, cosmic microwave background, large scale structure, indicate that the dark 
energy closely mimics a cosmological constant, and that the crossing of the "phantom 
line" , i.e. equation of state less than -1 is currently allowed. While these results represent 
a severe constraint for minimally-coupled quintessence models, a non-minimal coupling 
has been shown able to encompass it, by virtue of the enhanced dynamics imprinted at 
early time by the coupling itself. However, there is not a priori a particularly motivated 
form of this coupling, and the range of possibilities in which it can be selected is 
unlimited. From the cosmological point of view, however, it would be unphysical, 
as well as unpleasant, to rely on models suffering of the fine tuning problem for the 
initial configuration. Focusing on scalar-tensor theories as a viable model to connect 
Dark Energy and gravity, we have thus explored the possible choices of the coupling 
between a scalar field and the Ricci scalar, in order to select those ones giving rise to 
scaling solutions for the Klein Gordon equation, in the form (|12|). Our analysis resulted 
in selecting three classes of couplings, namely functions of the exponential form (122(1 . 
polynomial functions (|27)|) . and the exponential of polynomial (|2*7j) . depending on the 
coefficient m characterizing the scaling (jllj) and on the value of B on (J23J). 
Our analysis is complete, in the sense that it recovers all the possibilities to have a 
scaling behavior. Most importantly, it has been found that these solutions actually 
possess attractor properties, which is an extremely appealing feature in view of the old 
fine-tuning problem of minimally-coupled quintessence models. In particular, we have 
shown that all the scaling solutions of the form (jl2J) with n < 3 and 4 in the matter 
and radiation dominated era, respectively, represent attractors. Clearly, this enforces 
the case for extending the theory of gravity beyond general relativity, and opens a 
window on the solution of "initial values" problem. Cosmological models, characterized 
by the non-minimal couplings we have selected out in this papers, will deserve a further 
investigation. 
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